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Abstract 

We have performed numerical simulations of the 
interplanetary transport of solar cosmic rays. The 
particles form a coherent pulse within ~ 0.01 AU af- 
ter their injection. The gradual decrease of a pulse's 
speed and anisotropy can be understood in terms of 
an equilibrium between pitch-angle scattering and fo- 
cusing. The results should be useful for estimating 
times of particle injection. 

Introduction 

Much of the current research on solar cosmic rays 
aims to determine the time distribution of their emis- 
sion during solar flares. The major obstacle is that 
no spacecraft has approached closer than w 0.3 AU 
from the Sun, so interplanetary scattering has greatly 
broadened the observed temporal distributions. This 
problem is least severe when the detector is magneti- 
cally reasonably well-connected to the flare. If strong 
anisotropics are observed, they also aid the deter- 
mination of the scattering mean free path and the 
injection profile (Kallenrode, Wibberenz, and Hucke 
1992). 

The key phenomenon underlying the strong aniso- 
tropy observed for many well-connected events is that 
charged particles of a given energy tend to travel 
together as a highly anisotropic "coherent pulse," 
which can survive long enough to be observed near 
the Earth. Figure 1 shows stages in the formation 
of such a pulse very close to the Sun. This phe- 
nomenon was first systematically described by Earl 
(1974, 1976a, b) who used analytic approximations to 
describe their formation and decay. However, that de- 
scription was limited by the assumption of a constant 
focusing length, implying an exponentially decaying 
magnetic field strength. Here we re-examine the evo- 
lution of coherent pulses of solar cosmic rays using nu- 
merical simulations for a more realistic Archimedean 
spiral field (Parker 1958). 



For ions of ;> 20 MeV/n or electrons of ;> 10 keV, 
the transport along the interplanetary magnetic field 
is well described by the following equation (using the 
notation of Ng & Wong 1979): 
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where z is the arclength along the magnetic field, /i is 
the cosine of the pitch angle (v z /v), F is the distribu- 
tion function, defined as the number of particles per z 
per /i in a given magnetic flux tube, L is the focusing 
length, B / (dB /dz), and (p is the coefficient of pitch- 
angle scattering. The first term on the right hand side 
of eq. (1) expresses the streaming of particles along 
the magnetic field. The second term is due to pitch- 
angle scattering from irregularities in the magnetic 
field, and the final term is due to adiabatic focusing 
(Roelof 1969). 

We derive the focusing length, L(z), from an Archi- 
medean spiral field that makes a 45° angle with the 
radial direction at a radius of 1 AU, corresponding to 
a solar wind speed of « 400 km/s. The pitch-angle 
scattering coefficient, y(/i), is parameterized as 



^)=A| M r i (i- / i 2 ). 
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eq. (1) involves 
and not u or ( 



following Jokipii (1971). We express our results in 
terms of the scattering mean free path, A, which is 
given by 

A= (2- 9 )(4-g)T (3) 

Note that for a constant value of A 
only the distance traveled, s = vt. 
alone. We will therefore show results in terms of s 
instead of time. 

The numerical method used to solve eq. (1) was 
essentially the finite difference method of Ruffolo 
(1991), with some minor improvements. That study 
described tests of the code and successful fits to the 
observed intensity and pitch-angle distribution as a 
function of time for neutron-decay protons. More re- 
cently, this and two completely different numerical 
methods were shown to give very similar results (Earl 
et al. 1995). 
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Results 

Numerical simulations were performed for an ini- 
tial condition corresponding to an instantaneous in- 
jection of particles near the Sun at t = 0. For most of 
the simulations, we took the initial pitch-angle distri- 
bution to be concentrated at the highest /z-grid value, 
and used an absorbing boundary condition at z = 0, 
i.e., F(t > 0,/j, > 0,z = 0) was set to zero. (The 
injection site and inner boundary, z = 0, was set to 
r = 0.01 AU « 2Rq to avoid a divergence in the focus- 
ing term.) The outer boundary was set to the maxi- 
mum value of s — vt. Finally, the grid spacings were 
As = 0.02 AU, A/j, = 0.08, and Az = AfiAs = 0.0016 
AU except where noted. 

To examine the propagation of coherent pulses, we 
plot the average distance along the magnetic field of 
the particle distribution, (z), as a function of s = vt 
for selected values of A and q (Figure 2). The dashed 
line, (z) — s, shows where the particles would be if 
they all traveled directly along the magnetic field, 
with /i = 1 and v z = v. As expected, the progress 
along the magnetic field is fastest for the longest mean 
free path, corresponding to the weakest scattering. 

The rate of change of (z) with respect to s is 



ds 



(4) 



where 5 is the anisotropy averaged over the particle 
distribution, so {/i) is proportional to the pulse's prop- 
agation speed and anisotropy. Figure 3a shows the 
evolution of (fi) as a function of the average position, 
(z). The most noticeable trend is that (^) monotoni- 
cally decreases with position in all cases, as the pulse 
moves to regions where the focusing is weaker. This 
slowing in turn causes the curvature of the lines in 
Figure 2 away from the line of maximum speed. 

Another noticeable feature in Figure 3a is the con- 
vergence of curves for different q toward a universal 
curve for smaller (fi) values. This convergence occurs 
later for higher values of A and always appears at 
similar values of (/i), indicating that the convergence 
coincides with the decay of the coherent pulse toward 
a diffusive distribution. This can be verified by ex- 
amining a z = y/ ((z — {z)) 2 ) (Figure 3b). We note 
that the position at which a z ~ (z)/4 (long-dashed 
line) for each A corresponds roughly to the position 
of convergence, and to a value of (/i) « 1/3. We 
will therefore refer to (/i) ^ 1/3 as the coherent pulse 
regime and to (n) ^ 1/3 as the diffusive regime. 

The behavior of (/x) vs. (z) can be interpreted in 



terms of an equilibrium between the scattering and 
focusing terms in eq. (1) with the streaming term 
"turned off." [The scattering/focusing eigenfunc- 
tions defined by Earl (1976a) used a different form 
of eq. (1), and thus have a different physical inter- 
pretation.] Figure 4 shows simulation results for the 
equilibrium average of the pitch-angle cosine, {/J,) e q, 
and the equilibration distance, s eq , which is the dis- 
tance traveled over which equilibrium is achieved, or 
v times the longest relaxation time. The results for 
{/J,) eq are in good agreement with the analytic expres- 
sion 



eg 



(5) 

(fortunately, for the q values chosen here, these inte- 
grals can be solved in terms of simple functions). For 
large z (i.e., L ^> A, or the weak focusing limit), s eq 
approaches a g-dependent value (equal to or slightly 
larger than A) which agrees with the L — oo value of 
Bieber (1978). 

For very small values of z, where there is very 
strong focusing, s eq is also small, and the distribu- 
tion achieves pitch-angle equilibrium very quickly. At 
the same time, ((J,) e q ~ 1, so the distribution "re- 
laxes" rapidly (within ~ 0.01 AU) to form a coherent 
pulse with nearly maximal anisotropy and highly col- 
limated motion along the magnetic field (Figures 1 
and 5; a finer grid spacing was used for these simula- 
tions). Note that the spiky features in Figure 1, before 
pitch-angle equilibrium is reached, are due to the low 
scattering coefficient, ip (/i), near /j, = 1. Initial in- 
jections that were isotropic or highly focused yielded 
essentially the same distribution beyond z = 0.03 AU, 
with the former lagging by only 0.004 AU. A reflecting 
boundary condition gave virtually identical results. 

Farther from the Sun, we see that (n) e q rapidly de- 
clines with z (Figure 4). Since the equilibration dis- 
tance rapidly increases, the pitch-angle equilibrium is 
not maintained and (/z) declines less rapidly. For a 
narrow pulse, the deviation from pitch-angle equilib- 
rium evolves according to 



<W - (M)eg) 

ds 



dz 



(6) 



However, this equation gives a steady-state deviation 
which is insufficient to account for the difference be- 
tween (fi) and (n) eq in the diffusive regime. The re- 
maining deviation can be explained by considering 
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the width of the distribution. During the diffusive 
phase, the average of (/x) e9 over the distribution can 
deviate significantly from the value at (z). Thus we 
see that the convergence of (/x) as a function of (z) 
for different q's in the diffusive phase is due to 1) the 
convergence of (fj,) eq (z) toward a q- independent func- 
tion for which can be shown to be ^X/L(z), 
and 2) the similarity of a z as a function of (z) for the 
three values of q at each A (Figure 3b). 

Discussion 

We note that the transport equation used here ne- 
glects several well-known effects. Diffusion and drifts 
perpendicular to the field could reduce the intensity, 
especially for well-connected events, for which neigh- 
boring flux tubes have a lower intensity, but this 
should not significantly affect the propagation diag- 
nostics discussed here. Adiabatic deceleration and 
convection were included in numerical simulations by 
Ruffolo (1995) and were shown to noticeably affect the 
propagation for ions of ^<C 20 MeV/n and electrons of 
^ 10 keV. 

On the other hand, this approach does overcome 
the main limitation of the previous, analytic descrip- 
tions of coherent pulses (Earl 1976a, b), i.e., the as- 
sumption of a constant focusing length (B oc e~ z / L ). 
We find that several of their conclusions do not apply 
for the case of a realistic Archimedean spiral field. We 
find that for this case, the propagation speed is not 
constant, and does not exhibit distinct supercoherent 
and coherent regimes, but rather decreases steadily 
along with (fi) eq . Other results from our simulations 
show that the pulse width increases roughly linearly 
with time, in contrast with the result for an expo- 
nential field that a z oc i 1/2 (Earl 1974). Of course, 
these limitations do not diminish the importance of 
that seminal work. 

It is hoped that these results will be applicable to 
studies of the time profile of particle injection from 
solar flares, or comparisons of the time of injection 
for different particle species. When a strong coher- 
ent pulse is observed, one can estimate A from the 
anisotropy and the pulse width (Figure 3). A lack of 
consistency could indicate a finite injection width (if 
the pulse width is high) , or shock acceleration outside 
the corona (if the anisotropy is high) . One can derive 
the propagation time from A, given an assumption or 
measurement for q. Subtracting the propagation time 
from the time of maximum then yields an estimated 
time of injection. Note that conditions with A ~ 1 



AU, corresponding to the highest mean free paths re- 
ported in the inner heliosphcre, arc sometimes called 
"scatter-free" conditions, with the implication that 
particles freely stream along the field at their max- 
imum speed. However, our results show that even 
then, scattering significantly delays the arrival of the 
bulk of the pulse, while the "onset" is delayed some- 
what less (see also Kallcnrode and Wibberenz, 1990). 
In general, Figure 2 can be used to derive an accurate 
estimate of the propagation time. 
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Figure 2. Simulation results for the mean distance along the magne 
q — 1.0 (dotted lines), q — 1.5 (dashed lines), and q = 1.9 (solid lines) 
line corresponds to motion directly along the magnetic field. 

Figure 3. a) Mean pitch-angle cosine, (pi) = (v z )/v, and b) pulse wic 
(z). For an explanation of the curves, see Figure 1. Solid circles show pc 
that only curves for A = 1.0 AU are extended to s = 8.0 AU. 

Figure 4. a) Mean pitch-angle cosine, (//}, vs. mean distance, (z), up t 
for Figure 1, compared with {/J,) e q for A = 0.3 and q — 1.0 (o), 1.5 (®), a 
A and q as above. 

Figure 5. a) Mean pitch-angle cosine, (/i), vs. mean distance, (z), up t< 
(o) and highly focused (•) initial distributions, compared with {/i) eq (sc 



Figure 1. Three stages in the evolution of a coherent pulse: the distribution function, F (vertical direction), vs. distance 
along the magnetic field, z, and pitch-angle cosine, fi, for s = vt = 0.002, 0.01, and 0.02 AU, when A = 0.3 AU, q = 1.5, and 
F(t — 0) is concentrated at z = and uniform in /i. 



This preprint was prepared with AGU's IATfiJt macros v4. 
File 9506133 formatted February 1, 2008. 



r igure z 




r igure o 




r igure 4t 




(z) (AU) 



r igure o 



1.0 




0.05 



(z) (AU) 



